Generalized Schrieffer- Wolff Transformation of 2 Kondo Impurity Problem 
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We have carried out a generalized Schrieffer-Wolff transformation of an Anderson two-impurity 
Hamiltonian to study the low-energy spin interactions of the system. The 2 nd -order expansion 
yields the standard Kondo Hamiltonian for two impurities with additional scattering terms. At 
4" l -order, we get the well-known RKKY interaction. In addition, we also find an antiferromagnetic 
superexchange coupling and a correlated Kondo coupling between the two impurities. 

PACS numbers: 



I. INTRODUCTION 



The low-energy interactions of a system of magnetic 
impurities coupled to a sea of conduction electrons dis- 
play a rich variety of physical phenomena. One of the 
most "well-known many-body effects is screening of the 
local moment by the conduction electrons, first explained 
by Jun Kondo [lj ■ The Kondo effect was first observed in 
measurements of dilute magnetic impurities in metals , 
and is now ubiquitious in many systems including quan- 
tum dots (QD) 0, and magnetic atoms in surfaces 
probed by scanning tunneling microscopy 5]. 




FIG. 1: The magnetic impurities and the electron reservoirs 
are labeled 1 and 2, and L, M and R respectively. The tunnel- 
ing between the localized electron on impurity 1 and the con- 
duction electrons in the left and middle reservoirs are shown 
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In materials with a dense lattice of magnetic impuri- 
ties, there is an interesting interplay between the Kondo 
effect and magnetic ordering of the local moments, which 
is of particular interest in the heavy- fermion compounds. 
As the simplest manifestation of this competition, and 
as a first step towards understanding the lattice problem, 
considerable attention has been placed on the two Kondo 
impurity problem. The problem was first solved exactly 
in the particle-hole symmetric case using NRG jg and 
subsequently by the Bethe ansatz, conformal field theory 
, and bosonization Q . A second-order phase transition 
is found for the case of particle-hole symmetry, which typ- 
ically turns into a cross-over when potential scattering is 
introduced. 



A recent two quantum dot experiment by Craig et al. 
9] in the geometry shown in Fig. ^ showed a change 
in the Kondo resonance of one dot by tuning the sec- 
ond dot from an even to an odd occupation, ft has been 
suggested that this is due to RKKY coupling between 
the dots. We have looked at the effective inter-dot in- 
teractions obtained via a canonical transformation of the 
two impurity Anderson model. The purpose of the cal- 
culation is to extract the low-energy spin excitations of 
the system in order to predict the magnetic interactions. 
We expect such interactions to occur in STM systems of 
magnetic impurities on surfaces as well. 



II. ANDERSON MODEL OF QD 

The interactions of a quantum dot (QD) with its leads 
can be modelled using the Anderson Hamiltonian [Io| . 

|12|. For a two impurity system in the geometry as 
shown in Fig. 1, we propose a two impurity Anderson 
Hamiltonian. 
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Here, e a r is the kinetic energy of the conduction elec- 
trons in the leads (a = L, M and R), and t '% and i 2, - are 

n,k n,k 

the tunneling rates between the localized d-electrons at 
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sites 1 and 2, and conduction electrons at lead a. Ui and 
6j are the Coulomb charging energies and single-particle 
energies of the localized electrons in the two QDs. 

The low-energy spin interactions of the system are ob- 
tained from a generalized form of the Schrieffer- Wolff 
transformation, carried out on the Anderson two impu- 
rity Hamiltonian. The standard Schrieffer- Wolff trans- 
formation |13| of a a single Anderson impurity Hamilto- 
nian maps it to a single Kondo impurity problem using 
a canonical transform with an operator S. The general- 
ized Schrieffer- Wolff transformation, similar to a canoni- 
cal transformation for the Hubbard model carried out by 
MacDonald et al [l4|, is defined in the following manner. 

H = cxp 5 H cxp~ s 
S = ( 2 ) 

i 

where is of order O(j^r). Si is determined self- 
consistently at each level of the expansion as shown be- 
low. The usual Kondo Hamiltonian is obtained at 2 nd or- 
der in this expansion, and we carry out this expansion to 
A th order to obtain the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) and any other interactions at this order. Since 
we are only interested in the low-energy physics, we will 
discard terms that take the system out of the ground- 
state manifold, i.e. terms that cost energy U . This gives 



the following requirement at each level of expansion. 

(J -P G )H P G = Q (3) 

where P G is the Gutzwiller projector onto the ground 
state, which in this case is the singly-occupied d-state on 
both local sites. At first-order, the following condition is 
the defining equation for S^\ 

{I -Pa) (Hj + [S^,H ]j P G = (4) 

A straight-forward calculation then gives 




+ (tl«*y e** di a c aS a ) (5) 

Here, we have dropped the subscript n for the single- 
particle level in each QD as we are only interested in the 
highest level near Ep. is anti-hermitian as required 
so that H will be hermitian. The effective Hamiltonian 
at 2 nd -order is given by H = Hq + ifSW, Hi], where we 
find the Kondo terms for dot 1 and 2, shown in Eq. |S| 
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Here, Jl-'- and J2Z Z are the Kondo coupling at QD 

k,k' fc,fc' 

1 and 2 respectively. The next two terms represent po- 
tential scattering at the two sites, and are important at 
the quantum cross-over point between Kondo-dominated 
and RKKY-dominated behaviour. The expressions for 
J and K in terms of the original parameters are given 
in Appendix ^ The last two terms simply renormal- 
ize the single particle energies of the local d-electrons. 
There are also double occupancy and double vacancy 
terms that cost energy U, but are cancelled at 2 nd -order 
by [SW,Hq], as required by the Gutzwiller projection. 

At 3 rd -order, all the terms generated by the commu- 
tators of Hi n t, and take the system out of the 



ground state, and hence are required to cancel completely 
with [S^ 3 \Hq]. However, these terms generate terms at 
4 t,l -order that give rise to the well-known RKKY inter- 
action. We also find superexchange terms, and an in- 
teresting interaction that we characterize as a correlated 
Kondo effect. There are also additional contributions 
to these interactions that come from the commutators 
shown in Eq. IB2I We describe these interactions below. 
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III. ^"-ORDER TERMS 

The expansion was carried out order by order using 
a Mathcmatica program written to calculate the many 
commutators involved, including checks at each order. 
The first term we find at 4*' l -order is the well-known 
RKKY interaction between the two local sites. In ad- 
dition, we also find an anti-ferromagnetic superexchange 
term and a correlated Kondo term. The RKKY inter- 
action can be understood as a 2 nd -OTder perturbation 
on the Kondo ground state, with an electron-hole ex- 
cited state; hence the RKKY term is of 0(J 2 /E F ). The 
superexchange process arises from particle-particle and 
hole-hole excitations, resulting in an excited state that 
involves either double occupancy or double vacancy on 
either of the local sites. The electrons hop from both 
inpurity sites to an intermediate state in the conduc- 
tion sea, exchange, and hop back. This can only oc- 
cur when charge fluctuations are allowed at the impurity 
sites, hence the superexchange term is of 0(J 2 /U). In 
the strong Kondo limit (U — > oo), the superexchange in- 
teraction disappears, which is why superexchange is usu- 
ally not discussed for metallic Kondo systems. 

The correlated Kondo interaction arises from interme- 
diate states that involve an electron of the Kondo singlet 
hopping from one site to the other and back. Since it 
also involves doubly-occupied or doubly-vacant excited 
states on the impurity sites, it is also of 0(J 2 /U) and 
will disappear in the large-t/ limit. 

Making the assumption that the constants J, K and 
I are approximately constant over the energy range we 
are interested in, we make the approximation J™'" w J. 

k,k' 

I^'g is the interaction strength of the double-occupancy 
term and K~'g is the td renormalization energy. Looking 
at the form of J, we can approximate Is'g ~ an( i f° r 

J/2. The 

expressions for and K°t : - are also given in App. [5] 

k,k' k,k' 

The simplified expressions for RKKY and superex- 
change are then, 

j2 rfc F £) 4 

Hrkky = — — — [^0(^^12)^0(^/^12) 

Ef 107T 

+ Ji(k F Ri2)Y 1 (k F R 12 ))]Si-S 2 (7) 



the symmetric case where e<j U /2, K~ 'g 
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(k F L)4 f,h(k F R 12 ) 

U 



(2?r)2 



k F R\ 



Si-S 2 (8) 



Ji and Yi are respectively the first and second type Bessel 
functions of the z t,l -order, and R\ 2 is the distance be- 
tween the two impurities. 

The spatial dependence of RKKY and superexchange 
are plotted in Fig. |2 and both show the characteris- 
tic n/k F dependence. The ratio of Ef to U was cho- 
sen to be 1 in this plot. It must be understood that it 
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FIG. 2: Spatial dependence of RKKY and superexchange 
plotted with e d = -500 fieV, U = 4 meV, and E F = 36 
meV. RKKY is the red (long dashed) curve, superexchange 
is the green (short dashed) curve, and the blue (solid) curve 
is the sum of the two. The inset shows the experimentally 
relevant range of site separation. 



is only the sum of RKKY and superexchange that is a 
physically relevant interaction between the two Kondo 
impurities. Therefore, by changing the ratio of Ef to U, 
we can favour either the RKKY or superexchange con- 
tribution, and hence a more anti-ferromagnetic or more 
ferromagnetic interaction. 

In addition, we also found a correlated Kondo term at 
A th -order, which modulates the Kondo coupling on dot 1 
by varying the electron density on dot 2, and vice- versa. 
The expression for the term is, 
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The local impurity spins, S\ and S2, only couple to the 
conduction electrons in the middle reservoir, S m 

t ^^S-c r,i in this case because of the geometry 

shown in Figure 1. The intermediate excited states of 
the density-modulated Kondo process involve doubly- 
occupied and doubly- vacant states on the second dot. 
Hence, the coupling is of order 0(1/U), and depends on 
the density of the second dot. 

The spatial dependence of the correlated Kondo term 
is plotted in Fig. |3 and instead of a n/k F periodicity 
characteristic of RKKY, it has a 2ir/k F periodicity. 

In addition to the above terms, we also found density 
modulated RKKY and superexchange terms, and also 
additional scattering terms at 4 tft, -order. We will ignore 
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FIG. 3: Spatial dependence of correlated Kondo term plotted 
for €d = —500 /ieV, (7 = 4 meV, and Ef = 36 meV, which 
gives Jp = 0.22. The red (dashed) curve is the Kondo cou- 
pling, Jp, and the blue (solid) line is the Kondo coupling plus 
the correlated Kondo term. The inset, as in Fig. |5] shows the 
experimentally relevant range of site separation. Note that for 
this interaction, the zero crossing does not occur until nearly 
4/k F . 

the scattering terms as they are of order 0(J/Ep) weaker 
than the scattering terms obtained at 2™ d -order. The 
complete expressions in terms of the original parameters 
for RKKY and superexchange are given in App. IU1 



relevant spin interaction between the two impurities is 
actually the sum of both the RKKY and superexchange 
interactions. By changing the ratio of Ef to U and the 
distance between the two impurities, one could poten- 
tially tune the system through a phase transition from a 
Kondo singlet ground state to an impurity singlet ground 
state. 

Similar processes give rise to the correlated Kondo 
term that varies the Kondo coupling on dot 1, and hence 
changes the Kondo temperature, via modulation of the 
density on dot 2. This allows us to exert a non-local con- 
trol over the Kondo effect on the first impurity. Due to 
the short-range nature of this interaction, the two impu- 
rities have to be no more than 3 Xp apart for this effect 
to be observable. This can be achieved in STM experi- 
ments on Cu and CuN surfaces where cobalt atoms can 
be placed as near as 3.6 A apart. 

These effects can also be experimentally tested in 
quantum dot systems. In quantum dot systems, the 2 
DEG density can be modulated by a back gate voltage, 
thus changing Xp. This would allow us to observe the 
oscillation of the RKKY and superexchange interaction, 
and also to tune the exchange coupling between the two 
impurity spins. The occupation of the local state can 
also be varied using a separate gate voltage, V g , and the 
Kondo coupling can be measured using transport mea- 
surements through the first dot. 



IV. DISCUSSIONS 
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The low-energy spin interactions obtained from a gen- 
eralized Schrieffer- Wolff transformation on an Ander- 
son model of two magnetic impurities include not only 
the Kondo interaction at 2 nd -order and the well-known 
RKKY interaction at 4"*-order, but also an antiferro- 
magnetic superexchange interaction, which, unlike the 
RKKY interaction, cannot be obtained from a perturba- 
tion treatment of the Kondo interaction. The superex- 
change interaction arises from doubly-occupied and hole- 
hole excitations or double- vacant and particle-particle ex- 
citations; hence it cannot be obtained from the Kondo 
Hamiltonian that does not allow fluctuations in the oc- 
cupation of the localized states. Therefore, for systems 
where U is not too large compared to Ep, the physically 
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APPENDIX A: KONDO COUPLING 
DEFINITIONS 



The definitions for Jp'g ■> ^i'p 
below. 



and are given 
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APPENDIX B: HIGHER-ORDER S-OPERATORS 

The defining equation for was given in Eq. and 
the defining equations for the higher-order are, 



(I -Pa) ([SW,H ] + ^[SM,H I }))P G = (Bl) 
(T-P G ) (V 3 ),ff ] + ^« + = (B2) 

{I-Pg) {[S {a) -H q ] + \ \S (z) ,Hi} + 1 [S {2) , [S {2) ,H }] + , i [S™ , H x ]] + 



The detailed expressions for the higher-order Si- 
operators are not shown here as they are too tedious. 
The effective Hamiltonian up to 4 t,l -order is then given 
by the sum of all terms up to and including the 4 t,l -order 
expansion. 
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APPENDIX C: COMPLETE EXPRESSIONS 

The complete expression for superexchange in terms of 
the original parameters is given in Eq. IC1I 
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The complete expression for RKKY is shown in Eq. 

EM 



6 



Hrkky 



E 



.l,m/.l,m\* ,2,m/,2,m\, 



k kl ^ m '^ ^ ^"i,fcl 



fe v fel y fcl v fe ; 



E 



fe,fei 



e 2 e 



,fcl 



£ 1 \ e m,fc- e 2 



1 



£2 \ e 



+ 



,ki - £ 2 e m ,fe - e 2 - t/b e m>fc l - e 2 - t/ 2 ^ 
1 1 1 \ 



e ^ ^m.k ^ ^1 ^m.kl ^ t 



(C2) 



r 



The complete expression for the correlated Kondo term available from the authors upon request, 
is too involved and will not be reproduced here, but is 
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